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Abstract 

We study the time evolution of decaying particles in renormalizable 
models of Relativistic Quantum Field Theory. Significant differences 
between the latter and Non Relativistic Quantum Mechanics are found 
— in particular, the Zeno effect seems to be absent in such RQFT 
models. Conventional renormalization yields finite time behaviour in 
some cases but fails to produce finite survival probabilities in others. 

The analysis of the finite time evolution in relativistic quantum field the- 
ory (RQFT) and, in particular, the behaviour at short times of the survival 
probability of unstable particles in RQFT is a relevant matter as far as the 
foundations of quantum mechanics are concerned and has some peculiari- 
ties which make it to be worth studying. But before going to RQFT, we 
will briefly recapitulate the main features of the survival amplitude (SA) at 
short times in non relativistic quantum mechanics (QM henceforth); see, for 
instance, for a comprehensive review. 
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^E-mail: jolu@delta.ft.uam.es 



1 



Let P(t) (t > 0) the survival probability (SP) for an unstable system 
represented in QM by some normalized state | and let H be the Hamil- 
tonian governing the evolution of the quantum system. The SP is given by 
P(t) =| | exp(-itH) | | 2 (h = 1). Then, assuming the finiteness of the 
energy dispersion in the initial (unstable) state, i.e., 

(AE) 2 = (^\H 2 \^) - ((^\H\^)) 2 < oo, (1) 

as indeed happens in QM, one can show that a very short times 

P(t)~l-(AE)V (t->0+). (2) 

Notice that the last equation shows the violation of the exponential decay 
law at very short times. The quadratic behaviour in t gives rise to the so- 
called quantum Zeno effect, whose consequences and possible experimental 
detection have been profusely discussed (see || and references therein for 
some recent work on this effect). 

Another interesting point — which, however, is not the main subject of 
this paper — is the behaviour of the survival probability for large t values. 
One knows that if the hamiltonian spectrum is bounded from below, then the 
decay cannot be exponential, and P(t) > P exp (t), as t — > oo. Now, whereas 
that property (positivity of the energy) always holds in QM, the situation is 
less clear in RQFT. Let H ren be the hamiltonian of a certain model of RQFT, 
after having carried out standard ultraviolet renormalization. To the best of 
our knowledge, H ren has been rigorously constructed for superrenormalizable 
models, like, for instance, the A0 4 and Yukawa RQFT, in 1+1 dimensions, 
but presumably not for renormalizable models of the kind dealt with in this 
work. In the superrenormalizable RQFT models for which H ren exists, the 
latter is lower bounded [(|. In any case, it is not the lower boundedness of 
the hamiltonian but, rather, the finiteness of the energy dispersion in the 
initial (unstable) state what is relevant for the existence of the Zeno effect. 

Regarding the RQFT case, it has been pointed out || that the unrenor- 
malized survival amplitude of an unstable particle has, in any model of RQFT 
(superrenormalizable or renormalizable), a singular behaviour as t — > + 
and, in order to provide a cure for it, a new characteristic time has been 
introduced. Such singularities seem to be related to those already discussed 
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by Bogoliubov and Shirkov || for the Schrodinger state vector in RQFT. 
More recently, the present authors have shown that in a superrenormalizable 
model of RQFT (namely, a direct extension of the scalar 3 -theory to the de- 
cay of an unstable particle) the above short-time singularities of the survival 
amplitude can be fully eliminated through a conventional renormalization 
procedure 0. Here, we are going to analyze the time evolution of the SP 
of unstable particles in renormalizable models of RQFT. Now, one should 
be careful about vacuum polarizing interactions and ultraviolet divergent 
renormalization. Let us clearly state the problem 

We consider — with respect to a given reference frame — some "free" 
Hamiltonian Hq which describes particles for t < in RQFT in the Schrodinger 
picture. Also, let P the total three-momentum operator, with [P, Hq] = 0, 
and let | i,pi) the initial state representing the unstable particle at t — 0, with 
four- momentum = E{ = (mj +Pi 2 ) 1 ^ 2 , and renormalized mass 

Thus, \i,pl) is a common eigenstate of both Hq and P. If |0) represents the 
vacuum of H , then 



where a + (p~i) is the associated creation operator (in the Schrodinger picture). 
The interaction Hamiltonian Hj, which acts for t > and is responsible for 
the decay has, by assumption, the following properties: 



2. It gives rise to a renormalizable local RQFT. 

3. It involves a small coupling constant, in the sense that one can reliably 
apply perturbation theory. 

Upon going over to the interaction picture (ip) and after some straight- 
forward calculations, we get the basic formula enabling us to compute the 
survival amplitude (SA) in renormalized perturbation theory 




(3) 



1. It commutes with P. 




f 

Jo 




Jo 



dt n (0 1 Tia^Hjih) . . . Hj(t n )a + (0)] | 0) (4) 
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where, for simplicity, we have omitted ip in the operators, as well as pi in a 
and a + . Of course, T stands for time-ordered product. One can now evaluate 
A(t) by using (4) and Wick's theorem in the standard way. Then it will be 
useful to introduce the following factorized form of the SA 

A{t) = exp(-iEit)A(t), (5) 

where the "reduced" survival amplitude A(t) — > 1 for any t > if Hi — > 0. 
Now we have to deal with the divergences generically known as disconnected 
vacuum contributions (dvc). For the sake of brevity, we do not present here 
the details of the calculation (see 0] for a more elaborated treatment), which 
yields the following result 

A(t) = A dvc (t)A ph (t). (6) 

A p h{t), to be called hereafter the physical survival amplitude, is the sum of 
all perturbative contributions to A(t) in each of which a(t) is not contracted 
with a + (t = 0). Therefore all the contributions to A p h are free of infinite vol- 
ume divergences; all the latter are included in the factor A dvc (t). We should 
notice that A p h — > 1 as Hj — > 0. We shall concentrate on A p h. As already 
said, the above procedure has been applied to the case of a superrenormaliz- 
able (0 3 -like) interaction J7|, having shown that, after a standard ultraviolet 
mass renormalization, the survival amplitude is finite for all t and it behaves 
linearly in t as t — > (no Zeno effect, then). Nevertheless, the scope of such 
superrenormalizable interactions is rather limited, so that, here, we are going 
to deal with the case of renormalizable interactions, which is more relevant 
from a physical point of view. 

We shall consider a (unstable) relativistic spin 1/2 particle (i), with renor- 
malized mass m i; which can decay into two relativistic particles (a, b), with 
renormalized masses m a , rrib, so that: i) m; > m a + mt,, and ii) b is a spin-1/2 
fermion while a is a scalar spinless boson (we are considering all the particles 
neutral for simplicity; the general case is physically equivalent). Now H is 
the sum of three free Hamiltonians: one, corresponding to particle a being 
similar to those in the superrenormalizable case (kinetic energy of a scalar 
spin-0 boson) and the other two, of particles i and b, being the well-known 
for relativistic free Dirac fermions. For t > the interaction Hamiltonian in 
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the Schrodinger picture is 



Hi = Xf d 3 x [V b (x)Vi(x)Q a (x)] + h.c. - 5m t J Sx N[^i(x)%(x)] - 

(Z ia -I) j d 3 x [f i( 7 .9)*i - mi*i%] + o.c.t. (7) 

where A is a renormalized dimensionless coupling constant. The mass renor- 
malization counterterm for the unstable i-fermion is, to order A 2 

6m = (2^p J d ' k (k*-ml + ie)(k*-m 2 b +ie) (§) 

The term proportional to Z i2 — 1 is the standard contribution corresponding 
to the wave-function renormalization for particle i, and o.c.t. denote other 
renormalization counterterms which are not relevant here. The calculation 
now turns out to be much more elaborated than in the superrenormalizable 
case. One has, to order A 2 , 

A ph (t) - 1 = -=r u(pi, a)K ren u(pi, a) (9) 



K 



K ren = K + iSmit + i(Z i>2 - 1)(7-Pi - m i)t (10) 

= tttt^ / d%d A k 2 772 2 T'^TiT 6 2 i • J {3 H k i + h-pl)- 

(27r) & J (ki — m 2 a + ie)(k 2 — m b + ie) 

exp[it(Ei -k?i- kl)} - 1 exp[it(k°i + k Q 2 - Ei)] - 1 



i(Ei - ki - k° 2 ) i(k°i + k° 2 - Ei) 



(11) 



where e — > + and u, u denote normalized Dirac spinors for the unstable 

— * 

particle with polarization a. In K, we integrate over k 2 and, by adding — ie 
to Ei (which is a consistent procedure), we perform a residue integration over 
k 2 . These integrations give K as a piece linear in t plus another contribution, 
S. In S, we integrate over by residues, then obtaining 

-^ren = HTi ren + A 2 ^, (12) 



S ren - [dm, + (j.pi - mi)(Z it2 - 1)] = 
. A 2 f d%_ j.jpi -h)+m b | , AU 

V) 4 J kf-ml + ie' ( Pi - kiY — m\ + ie ~ lAbs ^en (IS) 
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d 3 h 1 [expitjEj- E a - E b )-l][E bl ° + (k 1 -pl).j + m b \ 
(2tt) 3 AE a E b 1 (E a + E b -Ei- ie) 2 

[exp it(Ej + E a + E b ) - 1] [-Erf + (fc~j - pi). 7 + rre 6 ] 
(£ a + E b + Ei - it) 2 



+ 



}(14) 



where E a = (m 2 a + ki) 1 / 2 , E b = \m\ + (pj - &!) 2 ] 1 / 2 . 



5 is the sum of two contributions each of which seems logarithmically 
ultraviolet divergent through naive power counting. Fortunately, H is fi- 
nite for any t > 0, since each divergence proportional to 7 cancels out by 
symmetric integration, while those proportional to 7 cancel out with one 
another. S ren is the renormalized self-energy of the unstable particle. The 
quantities A and AbsTj ren (the absorptive part of S ren ), defined uniquely 
through standard methods, are finite. We shall define the quantity Hi out 
of the first term on the right hand side of Eq. (14) as follows: i)we replace 
exp it(Ei — E a — E b ) — 1 by i sin t(Ei — E a — E b ), ii) if A = E a + E b — E^ we re- 
place [A-ieY 2 by -in5'(A). (Recall that (A-ie)~ 2 = PA~ 2 -in5'(A), where 
P is the principal part and 5' stands for the derivative of the 5 function.) This 
gives 5 = Hi + H 2 , where H 2 embodies all the remaining contributions to 5. 
Some direct computations yield the following properties: 

1) After cancelling out a linearly ultraviolet divergent piece proportional to 
7 (by symmetric integration), A bears the form Ai + 7 M A2 iA1 , Ai and A2 iM 
being logarithmically divergent quantities. 

2) By employing Srrii and the corresponding logarithmically ultraviolet diver- 
gent expression for Z ij2 — 1 to order A 2 (see, for instance, ||), all logarithmi- 
cally ultraviolet divergences cancel out in X = Srrii + ^.pi — rrii){Zi^ — 1) + A, 
which becomes finite. Moreover u(pi, a)Xu(pl, a) is real. 

3) Hi is linear and t, and, moreover, —tAbsY, ren + A 2 Hi = 

4) H 2 is ultraviolet finite 

We shall define H 3 out of H 2 as follows: in the expression yielding the 
latter, we replace exp it[Ei =F (E a + E b )] — 1 by cos t[Ei =f (E a + E b )] — 1. Then 
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we have H 2 = H 3 + H 4 , where 

1 r d 3 h coat(Ei - E a - E b ) - 1 o, ( C -w, 1 , 

(KT^ [-^7 + (h - Vd-l + m b ]} (15) 
while H 4 contains all the remaining contributions to S 2 . Then 

-f^ren = + A 2 ^4 + A 2 ^3 

The key reasons for the successive definitions of A, Ej, j = 1,2,3,4 are 
the above results 1)- 4) plus the following ones: 

5) u(jpi, a)[itX + A 2 H 4 ]w(]t*, a) is purely imaginary, 

6) u(pi, a)X 2 E 3 u(pi, a) is real 

Consequently, the physical survival probability up to order A 2 is 

777 ■ 

P(t) =| A ph (t) | 2 = 1 + 2\ 2 [u(p-a)~ 3 u(ti, o)\ V- (16) 

We shall now deal with a very important physical requirement, namely, 
the consistency between (16) and Fermi's golden rule. For that purpose we 
let t — > oo and recall that, then, g~ 2 [l — cos(tg)] — > nt5(q) for any real q. 
A glance at Eq. (15) shows that only the first fraction on its r.h.s. yields a 
non- vanishing contribution as t — > oo. Now, using well-known properties of 
Dirac matrices and spinors, one readily gets 

77? • 

P(t)~i--±rt (f-oo), (17) 

where T can be writen as 

r = 77^7 — / 7777^77777^7 ^ 2 + mi7Tlb 5 ^ ~ k ± ~ k 2)- 18 

167X2777,7 Eaik^Ebih) 

Notice that V is relativistically invariant and coincides exactly (as a straight- 
forward computation employing Dirac matrices and their traces show) with 
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the total decay rate of the unstable particle in its rest frame (pi = 0). By 
taking into account that rrii/Ei = (1 — v 2 ) 1 ! 2 , v} being the velocity of the 
decaying particle, and introducing its lifetime r = at rest, Eq. (17) can 
also be cast as P(t) — 1 — t(l — vfY^/r, which is simply showing the rela- 
tivistic fact that an unstable particle with velocity v~l in an inertial frame has 
a lifetime r(l — vl 2 )- 1 / 2 in that frame. Then, Eqs. (17-18) are consistent with 
the linearization of the exponential decay behaviour exp[— t(l — v 2 ) 1 ^ 2 ^ in 
an inertial frame; T, the decay rate in the rest frame, being given by Fermi's 
golden rule 

Next, we shall investigate the behaviour of the survival probability in the 
rest frame (pi = 0) for small positive values of t. For that purpose, we write 
S 3 = F + + 7°Y1 (noticing that the contribution proportional to 7 vanishes by 
symmetric integration). We are interested in the behaviours of Y± for small 
positive t. Such behaviours are determined by large values of | k± | . It suffices 
to neglect all masses in Y + , and one readily obtains Y + ~ — m^t / (IQti) . The 
treatment of Y"_ requires more care, due to an additional factor Eb(ki) which 
appears with opposite signs in both contributions in the r.h.s. of (15) and, 

— # 

hence, leads to cancellations at large | k± |. One finds Yl ~ — trrii/ (32%). 
Then, for very small values of t, the survival probability is given by 

A 2 

P(t) - 1 - — (2m b + mi)t (19) 

107T 

The linear behaviour for very small values of t appearing in (19) shows that 
the Zeno effect is not present. 

We shall now consider another renormalizable model, similar regarding 
physical description and notations to the one previously described, except 
for the following crucial differences: 

a) The unstable particle i is a spinless boson which decays into two spin 1/2 
(Dirac) fermions (a, b) (again, all particles are neutral for simplicity). 

b) For t > 0, the interaction hamiltonian is 

H T = X I d 3 (x) [* (x)tf 6 (x)$i(x)] + h.c. + c.t. (20) 
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The counterterms (c.t.) include both the mass and wave function renormal- 
ization counterterms for the decaying boson (see, for instance, [[J for their 
detailed expressions) plus another renormalization counterterms which are 
not relevant for our calculation. We shall limit ourselves to study the physi- 
cal survival amplitude A ph (t) for the decaying boson in its rest frame (pi = 0). 
A computation entirely similar to the one carried out for the previous model 
yields 

A ph (t) ~ 1 + ^{iXt + A 2 ~ 4 + A 2 2 3 ). (21) 

The quantities X and S 4 are the respective analogues of those for the previous 
model: both itX and S 4 are pure imaginary. X is the ultraviolet finite 
contribution which results after cancellations between: i) divergent terms 
arising from the actual analogue of K (Eq. (11) for the previous model) and 
the mass and wavefunction renormalization counterterms for the decaying 
boson, ii) finite contributions linear in t (analogues of —t AbsT, ren and A 2 5i 
in the previous model). We shall concentrate in the relevant piece, which is 
real and, hence, contributes to the survival probability P(t) = \A ph (t) | 2 : 

cost(E a + E b - rrij) - 1 cos t(E a + E b + rrii) - 1 

[ (E a + E b - mi y + (E a + E b + mi y u [ } 

where E afi = {m\ b + k x ) 1 ' 2 . 

The crucial difference with respect to the previous model is that, now, H 3 
is ultraviolet divergent for any finite value of t > 0, as a direct power counting 
shows. The same statement holds, in principle, for S 4 . To the best of our 
knowledge, and in spite of the fact that the present model is renormalizable, 
no cure seems to exist for such divergences. It is not hard to see that the 
source of the difficulties is related to the fact that the self-energy diagram 
for the boson is quadratically ultraviolet divergent. Therefore, it seems that 
if the relevant self-energy diagrams are, at most, logarithmically divergent, 
then one can obtain a finite survival probability for any t > 0. Contrarily, if 
higher order divergences are present, the survival probability does not exist, 
in the present context of strictly RQFT, for any finite value of t > 0, and, 
then, a different regularization procedure, such as the one introduced in || 
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— which, in a strict sense, lies outside the province of RQFT — should be 
employed in order to define the physical survival probability for finite times. 
In any case, one should note that it is formally possible to let t — > oo in (22), 
and employ [1 — cos(tg)]/g 2 — > irt5(q), as in the previous case. This formal 
limiting process in the divergent (real) quantity S 3 leads to a well defined 
limit for the survival probability P(t), namely, the actual analogue of Eq. 
(17) holds as well in the present model: the actual T coincides with the total 
decay rate of the unstable boson in its rest frame. 

Going back to the question of quadratic divergences and the non-existence 
then — in a strict RQFT context — of the survival probability at finite times, 
it is worth recalling that, in QED, the photon self-energy (vacuum polar- 
ization) diagram is quadratically divergent by naive power counting but, 
actually, it turns out to be only logarihmically divergent (before renormal- 
ization) because of gauge invariance (see [§], for instance) — a feature which 
bears some relationship to the behaviour of the and Z° self-energies in 
the Weinberg-Salam model. Therefore one could speculate on whether gauge 
invariance should play a crucial role to properly analyze time evolution, and 
in particular the survival probability at finite times of unstable particles in 
a RQFT context. 

As a final comment, we would like to emphasize that our analysis has 
relied on computations carried out in renormalized perturbation theory. On 
the other hand, several renormalizable models on RQFT (like the A0 4 one) 
are accepted to have a trivial continuum limit |1(J in a non-perturbative 



(renormalization group) framework. In the latter case, the peculiarities of 
survival amplitudes of unstable particles at finite times, as discussed in this 
work, would become irrelevant. Here again, (non-abelian) gauge invariance 
might come to the rescue, at least for QCD: it does make the renormalized 
theory non trivial in the continuum limit and one could speculate that, upon 
forcing all ultraviolet divergences to be logaritmic, time evolution would be 
well defined. 
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